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Abstract 

\Q . We show that Kruzhkov's theory of entropy solutions to multidimensional 

scalar conservation laws |Krj can be entirely recast in L 2 and fits into the 
general theory of maximal monotone operators in Hilbert spaces. Our ap- 
proach is based on a combination of level-set, kinetic and transport-collapse 
approximations, in the spirit of previous works by Giga, Miyakawa, Osher, 

> 

1 A short review of Kruzhkov's theory 

^ ■ First order systems of conservation laws read: 

'x. 

d t u + J2d*i(Qi(u)) = 0, 
i=i 

or, in short, using the nabla notation, 

d t u + V, • (Q(u)) = 0, (1) 

where u = u(t, x) G K m depends on t > 0, x £ K d , and • denotes the inner 
product in R d . The Qi (for i — 1,- • -,d) are given smooth functions from 
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R m into itself. The system is called hyperbolic when, for each r G R d and 
each U G R m , the m x m matrix J2i=i,d T iQi{U) can be put in diagonal form 
with real eigenvalues. There is no general theory to solve globally in time the 
initial value problem for such systems of PDEs. (See B DLI4 lDa| IMal ISe] for 



a general introduction to the field.) In general, smooth solutions are known 
to exist for short times but are expected to blow up in finite time. Therefore, 
it is usual to consider discontinuous weak solutions, satisfying additional 'en- 
tropy' conditions, to adress the initial value problem, but nothing is known, 
in general, about their existence . Some special situations are far better 
understood. First, for some special systems (enjoying 'linear degeneracy' or 
'null conditions'), smooth solutions may be global (shock free), at least for 
'small' initial data (see |Klj . for instance). Next, in one space dimension 
d — 1, for a large class of systems, existence and uniqueness of global weak 
entropy solutions have been (recently) proven for initial data of sufficiently 
small total variation jBBj . Still, in one space dimension, for a limited class 
of systems (typically for m — 2), existence of global weak entropy solutions 
have been obtained for large initial data by 'compensated compactness' ar- 
guments |Ta| IDi| ILPSj . Finally, there is a very comprehensive theory in the 
much simpler case of a single conservation laws, i.e. when m = 1. Then, 
equation (JTJ is called a 'scalar conservation law'. Kruzhkov |Krj showed that 
such a scalar conservation law has a unique 'entropy solution' u G L°° for 
each given initial condition uq G L°°. (If the derivative Q' is further assumed 
to be bounded, then we can substitute L\ oc for L°° in this statement.) An 
entropy (or Kruzhkov) solution is an L°° function that satisfies the following 
distributional inequality 

d t C(u)+V x -(Q c (u)) <0, (2) 

for all Lipschitz convex function C : R — > R, where the derivative of Q c is 
defined by (Q c )' = C'Q'. In addition, the initial condition u is prescribed 
in L} n „, namely: 

tajfK*,.) -«.(*)!* = 0, (3) 

for all compact subset B of R d . Beyond their existence and uniqueness, the 
Kruzhkov solutions enjoy many interesting properties. Each entropy solution 
u(t, ■), with initial condition « > continuously depends on t > in Lj oc and 
can be written T(t)uo, where (T(t),t > 0) is a family of order preserving 
operators: 

T(t)u > T(t)u , Vt > 0, (4) 
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whenever uq > uq. Since constants are trivial entropy solutions to (JT|), it 
follows that if uq takes its values in some fixed compact interval, so does 
u(t, ■) for all t > 0. Next, two solutions u and u, with uq — Uq G L 1 , are L 1 
stable with respect to their initial conditions: 



for all t > 0. As a consequence, the total variation TV(u(t, •)) of a Kruzhkov 
solution u at time t > cannot be larger than the total variation of its initial 
condition uo- This easily comes from the translation invariance of and 
from the following definition of the total variation of a function v. 



where || ■ || denotes the Euclidean norm on R d . The space L 1 plays a key 
role in Kruzhkov's theory. There is no LP stability with respect to initial 
conditions in any p > 1. Typically, for p > 1, the Sobolev norm | \u(t, •) | Iw 1 * 
of a Kruzhkov solution blows up in finite time. This fact has induced a 
great amount of pessimism about the possibility of a unified theory of global 
solutions for general multidimensional systems of hyperbolic conservation 
laws. Indeed, simple linear systems, such as the wave equation (written as a 
first order system) or the Maxwell equations, are not well posed in any LP 
but for p = 2 |Brnj . However, as shown in the present work, L? is a perfectly 
suitable space for entropy solutions to multidimensional scalar conservation 
laws, provided a different formulation is used, based on a combination of level- 
set, kinetic and transport-collapse approximations, in the spirit of previous 
works by Giga, Miyakawa, Osher, Tsai and the author [Brl. IBr2l IBr3l iBr4, 

iGMirmo] . 

2 Kruzhkov solutions revisited 

2.1 A maximal monotone operator in L 2 

Subsequently, we restrict ourself, for simplicity, to initial conditions uq{x) 
valued in [0, 1] and spatially periodic of period 1 in each direction. In other 
words, the variable x will be valued in the flat torus T d = R d /Z d . 

Let us now introduce: 
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1) the space L 2 ([0, 1] x T d ) of all square integrable functions 

(a, x) G [0, 1] x T d — >• Y(a, i) 61, 

2) the closed convex cone K of all Y e L 2 such that <9 a F > (in the sense 
of distributions), 

3) the subdifferential of K defined at each point Y G K by: 

dK{Y) = {Z e L 2 , J(Y - Y)Z dadx < , VY G K} , (7) 

4) the maximal monotone operator (MMO) (see jBrzj): 

Y^-q(a)-V x Y + dK(Y), (8) 

where q(a) = Q'(a), and the corresponding subdifferential equation |Brzj : 

G 9 t F + g(a) • V^y + ^y). (9) 

From maximal monotone operator theory |Brzj . we know that, for each 
initial condition Y G K, there is a unique solution Y(t,-) £ K to (jOJ, for 
alH > 0. More precisely, we will use the following definition (which includes 
the possibility of a left-hand side q G L 2 ([0, 1])): 

Definition 2.1 Y is a solution to 

q (a) G d t Y + q(a)-V x Y + dK(Y), (10) 

ura'£/i initial value Y G if and left-hand side q^ G £ 2 ([0, 1]), i/: 

1) t —> Y(t, ■) G L 2 zs continuous and valued in K , with Y(0, ■) = Y , 

2) Y satisfies, in the sense of distribution, 

— J \Y - Z\ 2 dadx < 2 J (Y - Z)(q (a) - d t Z - q(a) ■ V x Z)dadx, (11) 
for each smooth function Z(t, a, x) such that d a Z > 0. 

Proposition 2.2 For each Y G K, and q Q G L 2 ([0, 1]) ; there is a unique 
solution Y to MU\) in the sense of Definition \2. il If both Y and go belong to 
L°° , then we have for all t > 0: 

-tsup(-g )+ + mfy < Y(t,-) < sup Y + t sup(g )+. (12) 
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IfV x Yo belongs to L 2 , then so do dtY(t, •) and V x Y(t, ■) for all t > 0. Two 
solutions Y and Y to (EJ) (with different left-hand side go o,nd %) are L 2 
stable with respect to their initial conditions Y and Y in K: 

\\Y(t r )-Y(t,-)\\ L * < \\Y -Y \\ L 2+t\\q -q \\ L 2. (13) 

for all t > 0. This is also true for all p > 1, when both Y — Y and go _ go 
belong to L p : 

\\Y(t, ■) - Y(t, -)\\ LP <\\Y - Y \\ LP + t\\q - q \\ LP . (14) 

For the sake of completeness, a brief proof of these (standard) results will 
be provided at the end of the paper. 

2.2 The main result 

Our main result is 

Theorem 2.3 Let Y = Y(t, a, x) be a solution to the subdifferential equation 
(0|) with initial condition Yq G L°° , with d a YQ > 0. Then, 

u(t,y,x)= / H(y — Y(t, a, x))da, (15) 
Jo 

defines a one parameter family (parameterized byyER) of Kruzhkov solution 
to (QJ) ; valued in [0, 1]. In addition, all Kruzhkov solutions, with initial values 
in L°° , can be recovered this way (up to a trivial rescaling). 

Let us rapidly check the last statement of our main result. We must show 
that any Kruzhkov solution U(t,x) with initial condition Uq(x) valued in L°° 
can be recovered from a solution to (JHJ). To do that, according to the first 
part of the theorem, it is enough to find an L°° function Y (a,x) such that 
d a Y > and 

U (x) = / H(y - Y (a,x))da, 
Jo 

for some y £ K, say y — 1. This is always possible, up to rescaling, by 
assuming: 

r < U (x) < 1 — r 
for some constant r > 0. Indeed, we set 

Uo(y, x) = max(0, min(l, y U (x))) 
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so that Uq(x) — Uq(1, x) and d y u > 0. 

Then, for each fixed x, we solve Uo(y, x) = a by y = Y (a,x), setting: 

Y (a,x) = -^, VaG [0,1], Vx G T d , 
Uq{x) 

so that 

Uo(y,x)= H(y -Y (x,a))da. 
Jo 

(Notice that F is valued in [0, r -1 ].) Finally, according to the first part of 
the theorem, we get 

U(t,x)= f 1 H(l-Y(t,x,a))da, 
Jo 

where Y is the solution to (jHJ) with initial condition Y . 

2.2.1 Remark 

Notice that, for all t > 0, the level sets of Y and U are related by: 

{(a,x), U(t,x) > a} = {(a, x), Y(t,a,x) < 1}. 

Thus, the method of construction of Y out of Uq and the derivation of 

U(t,x) from Y(t, a, x) can be related to level-set methods in the spirit of 

|FSS[ IGilt IQFt ITGOj . This is why we may call 'level-set formulation' of 
scalar conservation law (J^J the subdifferential equation given by (JOJ) 

2.2.2 Remark 

The solutions (t,x) — ► u(t,y,x), parameterized by y G K, are automati- 
cally ordered in y. Indeed, d y u > immediately follows from representa- 
tion formula (|15|) . This is consistent with the order preserving property of 
Kruzhkov's theory (as explained in the first section). 

2.3 A second result 

The function u(t,y,x), given by (JT3J), can also be considered as a single 
Kruzhkov solution of a scalar conservation law in the enlarged (1 + d) dimen- 
sional space M x T d , namely 

d t u + d y (Q (u)) + V, ■ (Q(u)) = 0, (16) 
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with (y, x) G K x T d , provided: 

1) Q is zero, 

2) the initial condition u (y, x) is valued in [0, 1] and d y u > 0. 
Furthermore, it turns out that, if we add the left-hand side qo(a) = Q'o( a ) to 
©, so that we get ((TDJ) : 

qo(a) G d t Y + q(a)-V x Y + dK(Y), 

and solve for Y, then the corresponding u given by ()15|) is a Kruzhkov solution 
to flBjl. 

As a matter of fact, our proof will be done in this larger framework. We 
assume that q , q and Y are given in L°°, for simplicity. Without loss of 
generality, up to easy rescalings, we may assume that both q and Y are 
nonnegative, which simplifies some notations. 

Theorem 2.4 Assume that qo and q are given in L°° , with go > 0. Let 
Y = Y(t,a,x) be a solution to the subdifferential equation il(J\) . with initial 
condition Y G L°° , Y > and d a Y > 0. Then, 

u(t,y,x)= / H(y — Y(t,a,x))da, (17) 
J o 

is the unique Kruzhkov solution to Mb\) with initial condition: 

Uo(y,%)=[ H(y - Y (a, x))da. (18) 

JO 

In addition, Y is nonnegative and can be recovered from u as: 

/•CO 

Y(t,a,x)= / H(u(t,y, x) — a)dy. (19) 
Jo 

Before proving the theorem, let us observe that the recovery of Y from u 
through ()19j) is just a consequence of the following elementary lemma which 
generalizes (in a standard way) the inversion of a strictly increasing function 
of one real variable: 

Lemma 2.5 Let: a G [0,1] — > Z(a) G M+ with Z' > 0. We define the 
generalized inverse of Z : 

v(y) = [ H(y - Z(a))da, Wy G M. 
Jo 
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Then v' > 0, H(y - Z(a)) = H(v(y) - a) holds true a. e. in (a, y) G [0, 1] x K 

and: _ 

Z(a) = / H(a-v(y))dy. 
Jo 

In addtion, for a pair (Z,v), (Z,v) of such functions, we have the co-area 
formula: 

/ \Z(a)- Z(a)\da= / / \H(y - Z(a)) - H(y - Z(a))\dyda (20) 
Jo Jo Jo 

1 /"GO /"OO 

/ \H(v(y) - a) - H(v(y) - a)\dyda= \v(y) - v(y)\dy. 
o Jo Jo 

To recover (|T^|) . we notice first that d a Y > follows from the very def- 
inition E3] of a solution to (jlOj) . Next, Y" > follows from (|T2*|) and the 

assumptions q > 0, Yo > 0. Then, we apply lemma 12.51 for each fixed 
x G T d and t > 0, by setting Z(a) = Y(t, a, x) and u(t, y, x) = v(y). 



2.3.1 Remark 

The function f(t,a,y,x) = H(y — Y(t,a,x)) = H(u(t,y,x) — a) valued in 
{0, 1} is nothing but the solution of the Lions- Perthame-Tadmor |LPTj 'ki- 
netic formulation' of (fTBj) . which satisfies: 

dtf + qo(a)d y f + q(a) ■ V x f = <9 ay u, 

for some nonnegative measure n(t,a,y,x). 



2.3.2 Remark 

As already mentioned, the solutions of ([T0|) enjoys the LP stability property 
with respect to initial conditions (fTljl . not only for p = 2 but also for all 
p > 1. The case p = 1 is of particular interest. Let us consider two solutions 
Y and Y of ()10|) and the corresponding Kruzhkov solutions u and u given by 
Theorem 12.41 Using the co-area formula (J2*U|) . we find, for all t > 0, 

/ / \u(t,y,x) -u(t,y,x)\dxdy = 

JR JT d 




H(u(t, y, x) — a) — H(u(t, y, x) — a)\dadxdy 
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= / / / \H(y — Y(t, a, x)) — H(y — Y(t, a, x))\dadxdy 
Jo Jm Jf d 

I \Y(t, a, x) — Y(t, a, x)\dxda < / / \Y (a, x) — Y (a, x)\dxda 

JT d JO JT d 

\u (y,x) - u (y,x)\dxdy. 



Thus, Kruzhkov's L 1 stability property is nothing but a very incomplete 
output of the much stronger L p stability property provided by equation (fTU|) 
for all p > 1 . 

2.3.3 Remark 

As a matter of fact, in Theorem 12 A\ it is possible to translate the LP stabil- 
ity of the level set function Y in terms of the Kruzhkov solution u by using 
Monge-Kantorovich (MK) distances. Let us first recall that for two proba- 
bility measures \x and v compactly supported on R D , their p MK distance 
can be defined (see (53] for instance), for p > 1, by: 



SliHiv) = SU P / 4>(x)d/j(x) + J ifj(y)du(y) 1 

where the supremum is taken over all pair of continuous functions (f) and i[> 
such that: 

4>(x) + ip(y) < \x - y\ p , Vx,?/GK D . 
In dimension D = 1, this definition reduces to: 



6Jfj,,u) = \ \Y 



LP, 



where Y and Z are respectively the generalized inverse (in the sense of Lemma 
12. 5j) of u and v defined on R by: 

u (y) = A*([-oo, y}), v(y) = u([-oo, y]), My G R. 

Next, observe that, for each x G T d , the y derivative of the Kruzhkov solution 
u(t, y, x), as described in Theorem 12.41 can be seen as a probability measure 
compactly supported on R. (Indeed, d y u > 0, u = near y = —oo and u — 1 
near y = +oo.) Then, the L p stability property simply reads: 



5 p (dyu(t,-,x),d y u(t,-,x))dx < / 5 p (d y u (-,x),d y u (-,x))dx. 

We refer to [BBL and |CFLj for recent occurences of MK distances in the 
field of scalar conservation laws. 
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3 Proofs 

Let us now prove Theorem 12 . 41 ( which contains the first part of Theorem 12. 31 as 
the special case q Q = 0). The main idea is to provide, for both formulations 
(llfijl and the same time-discrete approximation scheme, namely the 

'transport-collapse' method |Brl| IBr2|. IBr3|. IGMj , and get the same limits. 

3.1 A time-discrete approximation 

We fix a time step h > and approximate Y(nh,a,x) by Y n (a,x), for each 
positive integer n. To get Y n from F„_i, we perform two steps, making the 
following induction assumptions: 

d a Y n ^>0, 0<F n _! <supr + (^-l)/isupg , (21) 
which are consistent with our assumptions on Y . 

Predictor step 

The first 'predictor' step amounts to solve the linear equation 

d t Y + q{a) ■ V X Y = q (a) (22) 

for nh — h < t < nh, with Y n -\ as initial condition at t — nh — h. We exactly 
get at time t = nh the predicted value: 

Y*(a, x) = Y n -i(a, x - h q(a)) + h qo(a). (23) 

Notice that, since q is supposed to be nonnegative, the induction assumption 
(12 1[) implies: 

< Y* < sup Y + nh sup q . (24) 

However, although d a Y n _i is nonnegative, the same may not be true for d a Y*. 
This is why, we need a correction step. 

Rearrangement step 

In the second step, we 'rearrange' Y* in increasing order with respect to 
a G [0, 1], for each fixed x, and get the corrected function Y n . Let us recall 
some elementary facts about rearrangements: 
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Lemma 3.1 Let: a G [0, 1] — > X(a) G R+ an L°° function. Then, there is 
unique L°° function Y : [0, 1] — > R + , such that Y' > and: 

f 1 H(y -Y(a))da= [ H(y-X(a))da, \/y G K. 
Jo Jo 

We say that Y is the rearrangement of X . In addition, for all Z G L°° such 
that Z' > 0, the following rearrangement inequality: 

J \Y(a) - Z(a)\ p da < J \X(a) - Z(a)\ p da. (25) 

holds true for all p > 1 . 

So, we define Y n (a, x) to be, for each fixed x, the rearrangement of Y*(a, x) 
in a G [0, 1]: 

d a Y n >0, f 1 H(y-Y n (a,x))da= [ H(y -Y*(a,x))da, Vt/GK. (26) 
Jo Jo 

Equivalently, we may define the auxiliary function: 

u n (y,x)= [ H(y-Y*(a,x))da, Vt/6 1, (27) 
Jo 



i.e. 



and set: 



u n (y,x)=[ H(y - h q (a) -Y n -i(a,x - h q(a)))da, (28) 
Jo 

poo 

Y n (a,x)= H(a-u n (y,x))dy. (29) 
Jo 

At this point, Y n is entirely determined by F„_i through formulae (j2HI), (J2SJ), 
or, equivalently, through formulae ()28p . ()29jl . Notice that, from the very defi- 
nition (pHjl of the rearrangement step, u n , defined by (J2ZJ), can be equivalently 
written: 

u n (y,x)= H(y -Y n (a,x))da. (30) 
Jo 

Also notice that, for all function Z(a,x) such that d a Z > 0, and all p > 1: 

y |F n (a, sc) — Z(a, x)| p <ia(ia; < / |F r *(a, s) — Z(a, a;)| p <ia(ix (31) 

11 



L 2 formulation of scalar conservation laws 



follows from the rearrangement inequality (|25jl. Finlly, we see that d a Y n > 
is automatically satisfied (this was the purpose of the rearrangement step) 
and 

< Y n < sup Y + nh sup go- 
follows form (|24j) (since the range of Y* is preserved by the rearrangement 
step). So, the induction assumption (j2Tj) is enforced at step n and the scheme 
is well defined. 

3.1.1 Remark 

Observe that, for any fixed function of y, is the (generalized) 

inverse of Y n (a, x), viewed as a function of a, in the sense of Lemma 12*31 Also 
notice that the level sets {(a, y); y > Y n (a, x)} and {(a, y); a < u n (y,x)} 
coincide. 

3.2 The transport-collapse scheme revisited 

The time-discrete scheme can be entirely recast in terms of u n (defined by 
(l30|l). Indeed, introducing 

ju n (a,y,x) = H(u n (y,x) - a), (32) 

we can rewrite ()28|). (J2T?|) in terms of u n and jw n only: 

u n (y,x)= [ ju n -i(y - h q (a),x - h q(a),a)da. (33) 
Jo 

We observe that, formulae ()32p33j) exactly define the 'transport-collapse' 
(TC) approximation to (fTH|) . or, equivalently, its 'kinetic' approximation, 
according to [Br]] [Br2l |Br3l IHMj . 

3.3 Convergence to the Kruzhkov solution 

We are now going to prove that, on one hand, Y n (a, x) converges to Y(t, a, x) 
as nh — > t, and, on the other hand, u n (y,x) converges to u(t,y,x), where 
Y and u are respectively the unique solution to sub differential equation (fTUj) 
with initial condition Y (a, x) and the unique Kruzhkov solution to with 
initial condition 

u {y,x)= H(y -Y Q (a,x))da. (34) 
Jo 
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From the convergence analysis of the TC method |Brl[ IBr2[ IBr3[ IGMj . we 

already know that, as nh — > t, 

\u n (y, x) — u(t, y, x)\dydx — > 0, 

where u is the unique Kruzhkov solution with initial value uq given by (|34j) . 
More precisely, if we extend the time discrete approximations u n (y,x) to all 
t G [0, T] by linear interpolation in time: 

hi \ i st — nh nh + h — t 
u (t,y,x) = u n+1 (y,x) — h u n (y, x) , (35) 

then u h - u converges to in the space C°([0, T], L l (R x T d )) as h — >• 0. 
Following (|19|) . it is now natural to introduce the level-set function Y defined 
by (jl9j) from the Kruzhkov solution: 

/•CO 

Y(t,a,x) — / H(a — u(t,y,x))dy. 
Jo 

(Notice that, at this point, we do not know that Y is a solution to the 
subdifferential formulation (fTUjl !) Let us interpolate the Y n by 

, r h, \ ^ r , .t — nh , ^nh + h — t , , 

F h (t,a,x) = y„ +1 (a,x)— ^— + F n (a,a;) , (36) 

for all t G [nh, nh + /i] and n > 0. By the co-area formula (ffilj) . we have 

J \Y(t, a, x) — Y n (a, x)\dadx = J \u(t, y, x) — u n (y, x)\dydx. 



Thus: 



sup \\Y(t,-)-Y h {t,-)\\ Li < sup \\u(t,-)-u h (t,-)\\ L i^0 
te[o,T] te[o,T] 



and we conclude that the approximate solution Y h must converge to Y in 
(7°([0,T],L 1 ([0, 1] x T d )) as /i -> 0. Notice that, since the Y h are uniformly 
bounded in the convergence also holds true in C°([0, T], L 2 ([0, 1] x T d )). 

We finally have to prove that Y is the solution to the subdifferential 
formulation with initial condition Yq. 
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3.4 Consistency of the transport-collapse scheme 

Let us check that the TC scheme is consistent with the subdifferential formu- 
lation (flT)|) in its semi-integral formulation For each smooth function 
Z(t, a, x) with d a Z > and p > 1, we have 



J \Y n+ i(a,x) — Z(nh + h,a,x)\ p dadx 

< J \Y* +1 (a, x) — Z{nh + h, a, x)\ p dadx 
(because of property ffTTTf) due to the rearrangement step (J2I3)) 

= J \Y n (a,x — hq(a)) + hq (a) — Z(nh + h,a,x)\ p dadx 

(by definition of the predictor step 

= J \Y n (a, x) + h qo(a) — Z(nh + h, a, x + h q(a))\ p dadx 

= [ \Y n - Z(nh, -)\ p dadx + h Y + o{h) 



(by Taylor expanding Z about (nh,a,x)). Since the approximate solution 
provided by the TC scheme has a unique limit Y, as shown in the previous 
section, this limit must satisfy: 



Y - Z\ p dadx <p (Y-Z)\Y- Z\ p - 2 (q (a) - d t Z - q(a) ■ V x Z)dadx, 



in the distributional sense in t. In particular, for p = 2, we exactly recover 
the semi-integral version (fTT|) of (fTU|) . We conclude that the approximate 
solutions generated by the TCM scheme do converge to the solutions of |T0|) 
in the sense of Definition 12.14 which completes the proof of Theorem 12.41 



where: 
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4 Viscous approximations 

A natural regularization for sub differential equation (fTUj) amounts to sub- 
stitute a barrier function for the convex cone K in L 2 ([0,1] x T d ) of all 
functions Y such that d a Y > 0. Typically, we introduce a convex function 
: K — >] — oo, +oo] such that 0(r) = +oo if r < 0, we define, for all Y G K, 

= / (f)(d a Y)dadx } (37) 

and set $(Y) = +oo if y does not belong to K. Typical examples are: 

0(r) = -log(r), 0(r) = rlog(r), 0(r) = -, Vr > 0. 

T 

Then, we considered the perturbed sub differential equation 

G d t Y + q(a)-V x Y-q (a)+ed$(Y), (38) 

for e > 0. The general theory of maximal monotone operators guarantees 
the convergence of the corresponding solutions to those of (fTUj) as e — ► 0. It 
is not difficult (at least formally) to identify the corresponding perturbation 
to scalar conservation (|16p. Indeed, assuming 0(r) to be smooth for r > 0, 
we get, for each smooth function Y such that d a Y > 0: 

d$(Y) = -d a {4>'{d a Y)). 

Thus, any smooth solution Y to ()38|). satisfying (9 a F > 0, solves the following 
parabolic equation: 

c^y + g(a) ■ - q (a) = ed a ((f>' \d a Y)) . (39) 

Introducing, the function u(t,y,x) implicitely defined by 

u(t, Y(t, a, x),x) = a, 

we get (by differentiating with respect to a, t and x): 

(d y u)(t,Y(t,a,x),x)d a Y(t,a,x) = 1, 

(d t u)(t, Y, x) + (9 y u)(t, F, x)d t Y = 0, 
(V x .n)(t, F, x) + (9,u)(t, F, x)V x .F = 0. 
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Multiplying by (d y u)(t, Y(t, a, x), x), we get: 

- d t u - q(u) ■ V x u - q (u)d y u = ed y {(j)'{— )). (40) 

CJyUL 

In particular, in the case 0(r) = — logr, we recognize a linear viscous ap- 
proximation to scalar conservation law |Td^ : 

d t u + q(u) ■ V x u + q {u)d y u = ed 2 y u, (41) 

with viscosity only in the y variable. 

4.0.1 Remark 

Of course, these statements are not rigourous since the parabolic equations 
we have considered are degenerate and their solutions may not be smooth. 

4.0.2 Remark 

In the case of our main result, Theorem 12. 31 we have go — and the variable y 
is just a dummy variable in (JTJ). Thus, the corresponding regularized version 

- d t u - q{u) ■ V x u = ed y {(j>'{^-)). (42) 

includes viscous effects not on the space variable x but rather on the 'param- 
eter' y G K. This unusual type of regularization has already been used and 
analyzed in the level-set framework developped by Giga for Hamilton- Jacobi 
equations |Gi2j . and by Giga, Giga, Osher, Tsai for scalar conservation laws 

[uuirmu] . 

5 Related equations 

A similar method can be applied to some special systems of conservation 
laws. A typical example (which was crucial for our understanding) is the 
'Born-Infeld-Chaplygin' system considered in |Br4j . and the related concept 
of 'order-preserving strings'. This system reads: 

d t (hv) + d y (hv 2 - hb 2 ) - d x (hb) = 0, (43) 

d t h + dy(hv) = 0, d t (hb) - d x (hv) = 0, 
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where h, b, v are real valued functions of time t and two space variables x, y. 
In |Br4j . this system is related to the following subdifferential system: 

G d t Y - d x W + dK(Y), d t W = d x Y, (44) 

where (Y, W) are real valued functions of (t, a, x) and K is the convex cone 
of all Y such that d a Y > 0. The (formal) correspondence between ()43|) and 
(j4^j) is obtained by setting: 

h(t, x, Y(t, x, a))d a Y{t, x, a) = 1, 

v (t, x, Y(t, x, a)) = d t Y(t, x, a), b(t, x, Y(t, x, a)) = d x Y(t, x, a). 

Unfortunately, this system is very special (its smooth solutions are easily inte- 
grable). In our opinion, it is very unlikely that L 2 formulations can be found 
for general hyperbolic conservation laws as easily as in the multidimensional 
scalar case. 

6 Appendix: proof of Proposition 12.21 

In the case when q and Y Q belong to L°° and are nonnegative, we already 
know, from the convergence of the TC scheme, that there is a solution Y to 
(fTU|) . with initial value Y , in the sense of definition 12.11 From (}2"Tj) . we also 
get for such solutions, when q > and Y > 0, 

< Y(t, ■) < supF + £supg , Vt > 0. 

By elementary rescalings, we can remove the assumptions that both Yq and 
go are nonnegative and get estimate (|T2|). 

Let us now examine some additional properties of the solutions to fT0|) ob- 
tained from the TC approximations. First, we observe that, in the TC 
scheme, 

1) the predictor step (a translation in the x variable by h q(a) plus an addi- 
tion of h qo(a)) is isometric in all LP spaces, 

2) the corrector step (an increasing rearrangement in the a variable) is non- 
expansive in all L p . 

Thus the scheme is non-expansive in all L p ([0, 1] x T d ). More precisely, for 
two different initial conditions Y and Yq, and two different data go and go, 
all in L°°, we get for the corresponding approximate solutions Y n and Y n : 

\\Y n - Y n \\ LP < ||y n _! - F n _i|| L P + h\\q - q \\ LP . (45) 
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This shows that (fTljl holds true for all solutions of (fTUj) generated by the TC 
scheme. 

Since the scheme is also invariant under translations in the x variable, we get 
the following a priori estimate: 

\\V x Y n \\ LP < \\V X Y \\ LP . (46) 

Finally, let us compare two solutions of the scheme Y n and Y n = Y n+ i ob- 
tained with initial condition Y = Yi. Using (jJ5J), we deduce: 

J \Y n+ i(a, x) — Y n (a, x)\ p dadx < J |Y"i(a, x) — Yq(cl, x)\ p dadx 

— J \Y{(a, x)—Y (a, x)\ p dadx = J \Y (a, x—h q(a))+h q (a)—Y (a, x)\ p dadx. 

So we get a second a priori estimate: 

Wn+x ~ Y n \\ L p < {\\qo\\up + \\q\\L°°\\V x Y \\ LP )h. (47) 

Thus the solutions Y to (fTUj) obtained from the TC scheme satisfy the a 
priori bounds: 

II V*Y(t, 0| |iP < HV^oIIlp, (48) 

\\d t Y(t,-)\\ LP < ||go||Lp + ||g||L-||V x .F || LP . (49) 

Notice that, at this level, we still do not know if solutions, in the sense of 
Definition 12.11 exist when Y e K and go £ -^ 2 ([0;1]) are n °t m L°° and 
we know nothing about their uniqueness. This can be easily addressed by 
standard functional analysis arguments. 



Existence for general data 

Let Y e K and go £ L 2 ([0, 1]). We can find two Cauchy sequences in L 2 , 
labelled by k G N, namely Yq G K and g^ G L 2 ([0, 1]), made of smooth 
functions, with limits Y and go respectively. Let us denote by Y k the cor- 
responding solutions, generated by the TC scheme. Because of their L 2 
stability, they satisfy: 

sup \\Y k (t,-)-Y k '(t,-)\\ L2 < \\Y k -Y k '\\ LP + T\\ q k -q k '\\ L2 . 

te[o,r] 

So, Y k is a Cauchy sequence in C°([0, T], L 2 ) of solutions of (fTU|) in the sense 
of Definition ^. 11 with a definite limit Y. Definition ^. H is clearly stable under 
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this convergence process. So, we conclude that Y satisfies the requirements 
of Definition 12 . II and is a solution with initial condition Yq and left-hand side 
go- Notice that, through our approximation process, we keep the a priori 
estimates (jUJ),©, for general data q e L 2 ([0, 1]). 

Uniqueness 

Let us consider a solution Y to (|l(Jj). with initial condition Y G K and 
left-hand side go £ -^ 2 ([0;1]); i n the sense of Definition 12.11 By definition 
Y(t, •) G K depends continuously of t G [0, T] in L 2 . From definition (fTTj) . 
using Z = as a test function, we see that: 

j t \\Y(t,-)\\ 2 L2 < 2 1 F(t,a,x)g (a) dadx < \\Y(t, -)\\ 2 L2 + ||g|| 2 2 , 

which implies that the L 2 norm Y(t, •) stays uniformly bounded on any finite 
interval [0,T]. Thus, T > being fixed, we can mollify Y and get, for each 
e g]0, 1] a smooth function Y e , valued in K, so that: 

sup \\Y(t,-)-Y e (t,-)\\ L2 <e. (50) 
te[o,T] 

Let us now consider an initial condition Zq such that V x Zq belongs to L 2 . 
We know that there exist a solution Z to (|1()|) . still in the sense of Definition 
12.11 obtained by TC approximation, for which both d t Z(t, •) and V x .Z(t, ■) 
stay uniformly bounded in L 2 for all t G [0,T]. This function Z has enough 
regularity to be used as a test function in (fTTj) when expressing that Y" is 
a solution in the sense of Definition 12.11 So, for each smooth nonnegative 
function 9(t), compactly supported in ]0,T[, we get from (fTTj) : 

j {0\t)\Y - Z\ 2 + 20(t)(Y - ^)(g (a) - 9 t Z - g(a) • V a .Z)}dadxdt > 0. 

Substituting Y e for Y", we have, thanks to estimate ([5(J|h 

y {0'(t)|Y e - Z| 2 + 20(t)(y e - Z)(g (a) - - g(a) • V x Z)}dadxdt > -Ce, 

where C is a constant depending on 9, Z, q and q only. Since Z is also a 
solution, using Y e as a test function, we get from formulation (fTTj) : 

/ {0'(t)|Z - Y e | 2 + 29(t)(Z - Y e )(g (a) - 9 4 Y e - g(a) ■ V^Yjjrfadxdt > 0. 
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Adding up these two inequalities, we deduce: 

J {26'(t)\Y e -Z\ 2 +26(t)(Y e -Z)(d t (Y € -Z)+q(a)-V x (Y e -Z))}dadxdt > -Ce. 
Integrating by part in t G [0, T] and x G T d , we simply get: 

J e\t)\Y t - Z\ 2 dadxdt > -Ce. 
Letting e — > 0, we deduce: 

4- [\Y- Z\ 2 dadx < 0. 
dtJ 

We conclude, at this point, that: 

\\Y(t,-)-Z(t,-)\\ L 2 < \\Y -Z \\ L 2, Vi G [0, T] 

This immediately implies the uniqueness of F. Indeed, any other solution Y 
with initial condition Fo must also satisfy: 

\\Y(t,-)-Z(t,-)\\ L 2 < \\Y -Z \\ L i. 

Thus, by the triangle inequality: 

\\Y(t,-) - Y(t,-)\\ L 2 <2\\Y - Z \\ L2 . 

Since Z Q G K is any function such that V X Z belongs to L 2 , we can make 
| \Y — Z \ \ L 2 arbitrarily small and conclude that Y — Y, which completes the 
proof of uniqueness. 
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